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Abstract 

We present new measurements of CP-violation parameters in B° — > KgKgKg and K^n ^ decays 
based on a sample of 275 x 10 6 BB pairs collected at the T(45) resonance with the Belle detector 
at the KEKB energy-asymmetric e + e _ collider. One neutral B meson is fully reconstructed in one 
of the specified decay channels, and the flavor of the accompanying B meson is identified from its 
decay products. CP-violation parameters are obtained from the asymmetries in the distributions 
of the proper-time intervals between the two B decays. We obtain 

s k%k%k% = +1.26 ± 0.68(stat) ±0.18(syst), 
Ak%k%k% = +0.54 ± 0.34(stat) ±0.08(syst), 

S K o^ = -0.58±°;^(stat)±0.11(syst), 

A K o n0j = +0.03 ± 0.34(stat) ±0.11 (syst). 

All results are preliminary. 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 
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I. INTRODUCTION 



In the standard model (SM), CP violation arises from an irreducible phase, the 
Kobayashi-Maskawa (KM) phase [1], in the weak-interaction quark- mixing matrix. In par- 
ticular, the SM predicts CP asymmetries in the time- dependent rates for B° and B° decays 
to a common CP eigenstate fcp [2]. In the decay chain T(4S*) — > B°B° — > fcpft&g, where 
one of the B mesons decays at time tap to a final state fcp and the other decays at time 
t tag to a final state f tag that distinguishes between B° and B°, the decay rate has a time 
dependence given by 

e -\At\/r B o f n 

V(At) = \l + q- Ssm(Am d At) + Acos(Am d At) \. (1) 

4r B o L L J J 

Here S and A are CP- violation parameters, r B o is the B° lifetime, Am d is the mass difference 
between the two B° mass eigenstates, At = tcp — t tag , and the 6-flavor charge q — +1 
(— 1) when the tagging B meson is a -B (-B ). To a good approximation, the SM predicts 
S = —£f sin 20!, where = +1(— 1) corresponds to CP-even (-odd) final states, and ^4 = 
for both b — > ccs and b — > sgg transitions. Recent measurements of time- dependent CP 
asymmetries in B° — > J/ipKg [3] and related decay modes, which are governed by the 
b — > ccs transition, by Belle [4, 5] and BaBar [6] already determine sin 20i rather precisely; 
the present world average value is sin 20i = +0.726±0.037 [7]. This serves as a firm reference 
point for the SM. 

The phenomena of CP violation in the flavor- changing b — > s transition are sensitive to 
physics at a very high-energy scale [8]. Theoretical studies indicate that large deviations from 
SM expectations are allowed for time- dependent CP asymmetries in B° meson decays [9]. 
Experimental investigations have recently been launched at the two B factories, each of 
which has produced more than 10 8 BB pairs. Belle's measurements of CP asymmetries using 
the decay modes B° -> <j>K%, <f)K° L , K+R-K% f Q {980)K° s , rj f K° s , uK° s , and K° s tc°, which are 
dominated by the b — > sqq transition, yielded a value that differs from the SM expectation 
by 2.4 standard deviations when all measurements are combined [10]. To elucidate the 
difference, it is essential to examine additional modes that may be sensitive to the same 
b — > s penguin amplitude. 

Recently it was pointed out that in decays of the type -B — > P°Q°X°, where P°, Q° 
and X° represent spin-0 neutral particles that are CP eigenstates, the final state is a CP 
eigenstate [11]. The B° — > KgKgKg decay, which is a = +1 state, is one of the most 
promising modes in this class of decays. The existence of this decay mode was first reported 
by the Belle collaboration [12], and recently confirmed with larger statistics by the BaBar 
collaboration [13]. Since there is no u quark in the final state, the decay is dominated by the 
b — > sss transition. In this report, we describe the first measurement of CP asymmetries 
in the B° -> K^K^K^ decay based on a 253 fb" 1 data sample that contains 275 x 10 6 BB 
pairs. 

We also measure time- dependent CP violation in the decay B° — > Kg-ir '-/, which is not 
a CP eigenstate but is sensitive to physics beyond the SM [14, 15]. Within the SM, the 
photon emitted from a B° (B°) meson is dominantly right-handed (left-handed). Therefore 
the polarization of the photon carries information on the original 6-flavor; the decay is thus 
almost flavor-specific. The SM predicts a small asymmetry S ~ — 2(m s /m;,)sin 2<f>i, where 
mj, (m s ) is the 6-quark (s-quark) mass [15]. Any significant deviation from this expectation 
would be a manifestation of physics beyond the SM. Belle's previous measurement of time- 
dependent CP asymmetries in the decay B° — > K* ^ (K*° — > Kgir ) [10] required that 
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the invariant mass of the K s and tt° (M K o n o) be between 0.8 and 1.0 GeV/c 2 to select the 
K*° — > K S 7T° decay [14] . Recently it was pointed out that in decays of the type B° — > P°Q 0/ ~f 
(e.g. P° = Kg and Q° = ir°), possible new physics effects on the mixing-induced CP 
violation do not depend on the resonant structure of the P°Q° system [15]. In this report, 
we describe a new measurement of CP asymmetries extending the invariant mass range to 
0.6 < M K o w o < 1.8 GeV/c 2 , which includes the K*° mass region, based on the new proposal. 

At the KEKB energy-asymmetric e + e" (3.5 on 8.0 GeV) collider [16], the T(4S f ) is pro- 
duced with a Lorentz boost of = 0.425 nearly along the electron beamline (z). Since 
the B° and B° mesons are approximately at rest in the T(AS) center-of-mass system (cms), 
At can be determined from the displacement in z between the fcp and / tag decay vertices: 
At ~ {z CP - z tag )/{(3 1C ) = Az/(/3 7 c). 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K\ mesons and to identify muons 
(KLM). The detector is described in detail elsewhere [17]. Two inner detector configurations 
were used. A 2.0 cm radius beampipe and a 3-layer silicon vertex detector (SVD-I) were 
used for a 140 fb^ 1 data sample (DS-I) containing 152 x 16 6 BB pairs, while a 1.5 cm radius 
beampipe, a 4-layer silicon detector (SVD-II) [18] and a small-cell inner drift chamber were 
used for an additional 113 fb _1 data sample (DS-II) that contains 123 x 10 6 BB pairs for a 
total of 275 x 10 6 BB pairs. 

II. EVENT SELECTION, FLAVOR TAGGING AND VERTEX RECONSTRUC- 
TION 

A. B° -h. K%K%K% 

We reconstruct the B° -> K S K S K S decay with a K^K^K^' or K+- K+' Kf final 
state, where a it + it~ (7r°7r°) state from a K s decay is denoted as Kg" (^s°)- We use these 
notations whenever appropriate. Although the Kg~ Kg~ K s ° state suffers from a larger 
background and a lower vertex reconstruction efficiency, we include it because its total 
branching fraction is larger than that of the Kg~ Kg~ Kg~ mode; with B(K S — > -r + tt") = 
2/3 and B(K° S -> 7r°7r°) = 1/3, we obtain B{3K° S -> Kp K^~ K^~) = (2/3) 3 = 8/27 while 
B{3K° S -> Kg—Kg—Kg ) = 3(2/3) 2 (l/3) = 12/27. We do not include final states with two 
or three K s — > 7r°7r° since their products of efficiencies and branching fractions are small. 

Pairs of oppositely charged tracks that have an invariant mass within 0.012 GeV/c 2 of the 
nominal K s mass are used to reconstruct K s — > ir + ir~ decays. The tt + tt~ vertex is required 
to be displaced from the IP by a minimum transverse distance of 0.22 cm for high momentum 
(> 1.5 GeV/c) candidates and 0.08 cm for those with momentum less than 1.5 GeV/c. The 
direction of the pion pair momentum must also agree with the direction defined by the IP 
and the vertex displacement within 0.03 rad for high-momentum candidates, and within 0.1 
rad for the remaining candidates. When we find two good Kg~ candidates that satisfy the 
criteria described above, we apply looser selection criteria for the third Kg~ candidate: (1) 
the mismatch in the z direction at the K s vertex point for the two n^ 1 tracks should be less 
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than 5 cm (1cm when both pions have associated SVD hits); (2) the angle in the r-0 plane 
between the K s momentum vector and the direction defined by the Kg and the IP should 
be less than 0.2 rad for high-momentum candidates, and less than 0.4 rad for the remaining 
candidates. 

Photons are identified as isolated ECL clusters that are not matched to any charged 
track. To select Kg — > 7r°7r° decays, we reconstruct tt° candidates from pairs of photons with 

> 0.05 GeV, where is the photon energy measured with the ECL. The reconstructed 
7T° candidate is required to satisfy 0.08 GeV/c 2 < M 77 < 0.15 GeV/c 2 and p c J s > 0.1 GeV/c, 
where M 77 and p^ are the invariant mass and the momentum in the cms, respectively. The 
large mass range is used to achieve a high reconstruction efficiency Candidate Kg — > 7r 7r 
decays are required to have invariant masses between 0.47 GeV/c 2 and 0.52 GeV/c 2 , where 
we perform a fit with constraints on the Kg vertex and the 7r° masses to improve the 
7r°7r° invariant mass resolution. We also require that the distance between the IP and the 
reconstructed K s decay vertex be between —10 cm and 100 cm, where the positive direction 
is defined by the K s momentum. The Kg* candidate is combined with two good Kg~ 
candidates to reconstruct B° — > Kg~Kg~K s ° decay, where we only use aforementioned 
good Kg~ candidates. 

For reconstructed B — > K S K S K S candidates, we identify B meson decays using the 
energy difference AE = E B ms — E^ m and the beam-energy constrained mass M bc = 

\Z(-^bcam) 2 — (PB ns ) 2 > where i?bcam * s the beam energy in the cms, and Eg ns and p§" s are the 
cms energy and momentum of the reconstructed B candidate, respectively. The B meson sig- 
nal region is defined as \AE\ < 0.10 GeV for B° -> K%~ K%- K%~ , -0.15 GeV < AE < 0.10 
GeV for B° -> K^'K^Kf, and 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 for both decays^ 

The dominant background to the B° — > KgKgKg decay comes from e + e~ — > uu, dd, ss, 
or cc continuum events. Since these tend to be jet-like, while the signal events tend to 
be spherical, we use a set of variables that characterize the event topology to distinguish 
between the two. We combine modified Fox- Wolfram moments [19] into a Fisher discriminant 
T. We also use the angle of the reconstructed B° candidate with respect to the beam 
direction in the cms (9b)- We combine JF and cosOb into a signal [background] likelihood 
variable, which is defined as £ S i g [bk g ] = £ S i g [bk g ] (^ r ) x £ S i g [bk g ] (cos 9b) ■ We impose requirements 
on the likelihood ratio 7Z s /b = Ai g / (£si g + £bk g ) to maximize the figure-of-merit (FoM) 
defined as / \J + iV bkg , where (N hkg ) represents the expected number of signal 
(background) events in the signal region. We estimate N^p using Monte Carlo (MC) events, 
while iVbk g is determined from events outside the signal region. The requirement for 7Z s /b 
depends both on the decay mode and on the flavor-tagging quality, r, which is described in 
Sec. IIC. The threshold values range from 0.2 (used for r > 0.875) to 0.5 (used for r < 0.25) 
for the case with 3 good Kg~ candidates, and from 0.3 to 0.9 for other cases. 

We reject K S K S K S candidates that are consistent with B° -> XcoK° s -> (K S K%)K S or 
B° -> D°K° S -> (K S K%)K S . We keep candidate B° -> / (980)Arg -> (K S K° S )K S decays as 
they are also dominated by the b — > s transition. 

We use events outside the signal region as well as a large MC sample to study the 
background components. The dominant background is from continuum. The contributions 
from BB events are small. The contamination of B° — > XcoKg events in the B° — > K S K S K S 
sample is small (2.6%). The influence of this background is treated as a source of systematic 
uncertainty. Backgrounds from the decay B° — > D°K S are found to be negligible. 

Figure 1(a) [(c)] shows the M bc [AE] distribution for the reconstructed B° — > K S K S K S 
candidates within the AE [M bc ] signal regions after flavor tagging and vertex reconstruction. 
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The signal yield is determined from an unbinned two-dimensional maximum-likelihood fit 
to the AE-Mbc distribution. The fit region for the B° — > K$~ K$~ K$~ decay is defined as 
5.20 GeV/c 2 < M bc < 5.30 GeV/c 2 and -0.30 GeV < AE < 0.50 GeV, excluding the region 
5.26 GeV/c 2 < M bc < 5.30 GeV/c 2 and -0.30 GeV < AE < -0.12 GeV to reduce the effect 
of background from B decays. The fit region for the -B — > Kg^Kg^K^P decay is defined as 
5.22 GeV/c 2 < M bc < 5.30 GeV/c 2 and -0.40 GeV < AE < 0.40 GeV, excluding the region 
5.25 GeV/c 2 < M bc < 5.30 GeV/c 2 and -0.4 GeV < AE < -0.17 GeV. The K+~ K+~ K+~ 
signal distribution is modeled with a Gaussian function (a sum of two Gaussian functions) for 
M bc (AE). For the B° — > Kg~K^~Kg° decay, the signal is modeled with a two-dimensional 
smoothed histogram obtained from MC events. For the continuum background, we use the 
ARGUS parameterization [20] for M bc and a linear function for AE. The fits after flavor 
tagging yield 72 ±10 5°^ K$~ K£~ K£~ events and 16 ± 8 B° -> Kg~Kg~Kg° events for 
a total of 88 ± 13 B° — > Kg Kg Kg events in the signal region, where the errors are statistical 
only. 

B. B° -» K°7T°7 

Candidate K s — > 7r + 7r~ decays are selected with the same criteria as those used to select 
good Kg~ candidates for the -B — > K S K S K S decay, except that we impose a more stringent 
invariant mass requirement; only pairs of oppositely charged pions that have an invariant 
mass within 0.006 GeV/c 2 of the nominal K s mass are used. Candidate 7r° mesons are 
required to satisfy 0.118 GeV/c 2 < M 77 < 0.15 GeV/c 2 and p c J s > 0.3 GeV/c. The K° s ir° 
invariant mass, M^o^-o, is required to be between 0.6 and 1.8 GeV/c 2 . 

For prompt photons from the B° — > K s 7r 0, ~f decay, we select the photon with the largest 
E^ ms among photon candidates in an event and require 1.4 GeV < E^ ms < 3.4 GeV, where 
E^ ms is the photon energy in the cms. For the selected photon, we also require E 9 /E 2 5 > 0.95, 
where E9/E25 is the ratio of energies summed in 3 x 3 and 5x5 arrays of CsI(Tl) crystals 
surrounding the crystal at the center of the shower. Photons for candidate ir° — > 77 or 
V ~^ 77 decays are not used; we reject photon pairs that satisfy C n o > 0.18 or C v > 0.18, 
where C n o^ is a tt° (77) likelihood described in detail elsewhere [21]. The polar angle of the 
photon direction in the laboratory frame is required to be between 33° and 128° for DS-I, 
while no requirement is imposed for DS-II as the material within the acceptance of the ECL 
is much reduced for this dataset. 

Candidate B + — > K s n + ^ decays are also selected using a similar procedure to reconstruct 
the decay B° — > K^n ^. Candidate B + — > K° s tt + ^ and B° — > Kf^n ^ decays are selected 
simultaneously; we allow only one candidate for each event. The best candidate selection 
is based on the event likelihood ratio lZ s /b that is obtained by combining JF, which uses 
the extended modified Fox- Wolfram moments [22] as discriminating variables, with cos9 H 
defined as the angle between the B° meson momentum and the daughter Kg momentum in 
the rest frame of the K^ir system. We select the candidate with the largest lZ s /b- 

The signal region for the B° — > ^7r°7 decay is defined as —0.2 GeV < AE < 0.1 
GeV, 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 . We require TZ s/h > 0.5 to reduce the continuum 
background. 

The selection criteria described above are the same as those used for the previous time- 
dependent CP asymmetry measurement in the B° — > K* 0/ ~f (K*° — > K^ ) decay [10], 
except for the wider M K o n o range. 
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We use events outside the signal region as well as a large MC sample to study the back- 
ground components. The dominant background is from continuum. Background contribu- 
tions from B decays are significantly smaller than those from continuum, and are dominated 
by cross-feed from other radiative B decays including B + — > K s n + 'y, and charmless B de- 
cays. Background from other BB decays is found to be negligible. 

Figure 1(b) [(d)] shows the M bc [AE] distribution for the reconstructed B° — > Kg7v°j 
candidates within the AE [M bc ] signal region after flavor tagging and vertex reconstruction. 
The signal yield is determined from an unbinned two-dimensional maximum-likelihood fit 
to the AE-M hc distribution in the fit region defined as 5.20 GeV/c 2 < M bc < 5.29 GeV/c 2 
and —0.4 GeV < AE < 0.5 GeV. The B° — > Kgir ^ signal distribution is represented 
by a smoothed histogram obtained from MC simulation that accounts for the correlation 
between M bc and AE. The background from B decays is also modeled with a smoothed 
histogram obtained from MC events; its normalization is a free parameter in the fit. For the 
continuum background, we use the ARGUS parameterization for M bc and a second-order 
Chebyshev function for AE. The fit yields 105 ± 14 B° — > Kgir 0/ -f events, where the error is 
statistical only. For reference, we also measure the signal before vertex reconstruction and 
flavor tagging, and obtain 221 ± 21 events. 

C. Flavor Tagging 

The 6-flavor of the accompanying B meson is identified from inclusive properties of par- 
ticles that are not associated with the reconstructed B° fcp decay [23]. We use the 
same procedure that is used for the sin20i measurement [5]. The algorithm for flavor tag- 
ging is described in detail elsewhere [24]. We use two parameters, q and r, to represent 
the tagging information. The first, q, is already defined in Eq. (1). The parameter r is an 
event-by-event, MC-determined flavor-tagging dilution factor that ranges from r = for no 
flavor discrimination to r = 1 for unambiguous flavor assignment. It is used only to sort 
data into six r intervals. The wrong tag fractions for the six r intervals, Wi (I = 1,6), and 
differences between B° and B° decays, Awi, are determined from the data; we use the same 
values that were used for the sin 20! measurement [5] for DS-I. Wrong tag fractions for DS-II 
are separately obtained with the same procedure; we find that the values for DS-II, which 
are listed in Table I, are slightly smaller than those for DS-I. The total effective tagging 
efficiency for DS-II is determined to be e e g = X)f=i e «(l — 2w;) 2 = 0.30 ± 0.01, where q is 
the event fraction for each r interval determined from the J/ipK^ data and listed in Table I. 
The error includes both statistical and systematic uncertainties. 

D. Vertex Reconstruction 

The vertex position for B° — > KgKgKg and B° — > Kg7i 0, y decays is reconstructed using 
charged pions from K° s decays. A constraint on the IP is also used with the selected tracks; 
the IP profile is convolved with finite B flight length in the plane perpendicular to the z axis. 
Both charged pions from the K$ decay are required to have enough SVD hits to reconstruct 
a Kg trajectory with a sufficient resolution. The reconstruction efficiency depends both on 
the Kg momentum and on the SVD geometry. Efficiencies with SVD-II are higher than 
those with SVD-I because of the larger outer radius and the additional layer. 
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The /tag vertex determination with SVD-I remains unchanged from the previous pub- 
lication [10], and is described in detail elsewhere [25]; to minimize the effect of long-lived 
particles, secondary vertices from charmed hadrons and a small fraction of poorly recon- 
structed tracks, we adopt an iterative procedure in which the track that gives the largest 
contribution to the vertex x 2 is removed at each step until a good x 2 is obtained. 

For SVD-II [10], we find that the same vertex reconstruction algorithm results in a larger 
outlier fraction when only one track remains after the iteration procedure. Therefore, in this 
case, we repeat the iteration procedure with a more stringent requirement on the SVD-II hit 
pattern. The resulting outlier fraction is comparable to that for SVD-I, while the inefficiency 
caused by this change is small (2.5%). 

E. Summary of Signal Yields 

The signal yields for B° — > fcp decays, N sig , after flavor tagging are summarized in 
Table II. For the B° — ► K S K S K S decay, results both before and after vertex reconstruction 
are listed. The result for the B° — > Kgir ^ is obtained after vertex reconstruction. The 
signal purities are also listed in the table. 

III. RESULTS OF CP ASYMMETRY MEASUREMENTS 

We determine S and A for each mode by performing an unbinned maximum-likelihood 
fit to the observed At distribution. The probability density function (PDF) expected for 
the signal distribution, P sig (At; S, A, q, Wi, Awi), is given by Eq. (1) incorporating the effect 
of incorrect flavor assignment. The distribution is convolved with the proper-time interval 
resolution function i? sig (At), which takes into account the finite vertex resolution. 

The At resolution function for the B° — > K s n 0/ y decay is the same as that used in the 
previous analysis on the decay B° — > K* ^ (K*° — > K s ir°) [10]. It is based on the resolution 
function obtained from flavor-specific B decays governed by semileptonic or hadronic b — > c 
transitions, with additional parameters that rescale vertex errors. The rescaling parameters 
depend on the detector configuration (SVD-I or SVD-II), SVD hit patterns of charged pions 
from the K s decay, and K s decay vertex position in the plane perpendicular to the beam 
axis. These parameters are determined from a fit to the At distribution of B° — > J/ipK s 
data. Here the K s and the IP constraint are used for the vertex reconstruction, the B° 
lifetime is fixed at the world average value, and 6-flavor tagging information is not used so 
that the expected PDF is an exponential function convolved with the resolution function. 

We check the resulting resolution function by also reconstructing the vertex with leptons 
from J/vf) decays and the IP constraint. We find that the distribution of the distance 
between the vertex positions obtained with the two methods is well represented by the 
obtained resolution function convolved with the well-known resolution for the J/ip vertex. 
Finally, we also perform a fit to the B° — > J/tpK s sample with 6-flavor information and 
obtain Sj/^ K o = +0.68 ± O.lO(stat) and Aj/^, K o = +0.02 ± 0.04(stat), which are in good 
agreement with the world average values. Thus, we conclude that the vertex resolution for 
the -B — > KgTi ^ decay is well understood. 

The resolution function for the B° — > K S K° S K S decay is based on the same resolution 
function parameterization. The rescaling factor depends on the detector configuration (SVD- 
I or SVD-II) , SVD hit patterns of charged pions from the K s decay, K s decay vertex position 
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in the plane perpendicular to the beam axis, and the number of Kg~ decays used for the 
vertex reconstruction. 

We determine the following likelihood value for each event: 

P l = (Wei) J[f sig V sig (At')R sig (At t - At') 

+ (1 - /sig)nk g (A0i4k g (A^ - A*')]d(Af) 

+ f ol P ol (AU) (2) 

where P \(At) is a broad Gaussian function that represents an outlier component with a small 
fraction / Q i. The signal probability / s i g depends on the r region and is calculated on an event- 
by-event basis as a function of AE and M^ c for each mode. A PDF for background events, 
'Pbkg(At), is modeled as a sum of exponential and prompt components, and is convolved 
with a sum of two Gaussians i?bk g - All parameters in Pbkg(At) and i?bk g are determined by 
the fit to the At distribution of a background-enhanced control sample; i.e. events outside 
of the AE-M hc signal region. We fix r B a and Am d at their world-average values [26]. In 
order to reduce the statistical error on A, we include events without vertex information in 
the analysis of B° — > KgKgK^. The likelihood value in this case is obtained by integrating 
Eq. (2) over At;. 

The only free parameters in the final fit are S and A, which are determined by maximizing 
the likelihood function L = flj Pi{Ati\ S, A) where the product is over all events. 

Table III summarizes the fit results of S and A. We define the raw asymmetry in each 
At bin by (N q=+ i — N q= -i) / (N q=+ i + N q= -i), where N q=+1 ^_^ is the number of observed 
candidates with q — +1(— 1). Figures 2(a) and (b) show the raw asymmetries for the B° — > 
KgKgKg and KgTT 0/ ~f decays, respectively, in two regions of the flavor-tagging parameter r. 
While the numbers of events in the two regions are similar, the effective tagging efficiency 
is much larger and the background dilution is smaller in the region 0.5 < r < 1.0. Note 
that these projections onto the At axis do not take into account event-by-event information 
(such as the signal fraction, the wrong tag fraction and the vertex resolution), which is used 
in the unbinned maximum-likelihood fit. 

Table IV lists the systematic errors on S and A. The total systematic errors are obtained 
by adding each contribution in quadrature, and are much smaller than the statistical errors 
for all modes. 

To determine the systematic error that arises from uncertainties in the vertex recon- 
struction, the track and vertex selection criteria are varied to search for possible systematic 
biases. Small biases in the Az measurement are observed in e + e~ — > fi + fi~ and other con- 
trol samples. Systematic errors are estimated by applying special correction functions to 
account for the observed biases, repeating the fit, and comparing the obtained values with 
the nominal results. The systematic error due to the IP constraint in the vertex reconstruc- 
tion is estimated by varying (±10/xm) the smearing used to account for the B flight length. 
Systematic errors due to imperfect SVD alignment are determined from MC samples that 
have artificial mis-alignment effects to reproduce impact-parameter resolutions observed in 
data. 

Systematic errors due to uncertainties in the wrong tag fractions are studied by varying 
the wrong tag fraction individually for each r region. Systematic errors due to uncertainties 
in the resolution function are also estimated by varying each resolution parameter obtained 
from data (MC) by ±ler (±2<r), repeating the fit and adding each variation in quadrature. 
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Each physics parameter such as r B a and Am<2 is also varied by its error. A possible fit bias 
is examined by fitting a large number of MC events. 

Systematic errors from uncertainties in the background fractions and in the background 
At shape are estimated by varying each background parameter obtained from data (MC) 
by ±1(7 (±2<x). Uncertainties in the background B decay model are also considered for the 
B° — > KgTT 0/ ~f mode; we compare different theoretical models for radiative B decays and 
take the largest variation as the systematic error. 

Additional sources of systematic errors are considered for B decay backgrounds that are 
neglected in the PDF. We consider uncertainties both in their fractions and CP asymmetries; 
for modes that have non- vanishing CP asymmetries, we conservatively vary the CP- violation 
parameters within the physical region and take the largest variation as the systematic error. 
The effect of backgrounds from XcoKg in the B° — > KgKgKg sample is considered. 

Finally, we investigate the effects of interference between CKM-favored and CKM- 
suppressed B — > D transitions in the / tag final state [27]. A small correction to the PDF for 
the signal distribution arises from the interference. We estimate the size of the correction 
using the B° — > D*~£ + v sample. We then generate MC pseudo-experiments and make an 
ensemble test to obtain systematic biases in S and A. 

Various crosschecks of the measurement are performed. We reconstruct charged B meson 
decays that are the counterparts of the B° —> fcp decays and apply the same fit procedure. 
All results for the S term are consistent with no CP asymmetry, as expected. Lifetime 
measurements are also performed for the fcp modes and the corresponding charged B 
decay modes. The fits yield r B o and t b + values consistent with the world-average values. 
MC pseudo-experiments are generated for each decay mode to perform ensemble tests. We 
find that the statistical errors obtained in our measurements are all consistent with the 
expectations from the ensemble tests. 

A fit to the B° -> ^7T°7 subsample with 0.6 < M K ^ < 0.8 GeV/c 2 or 1.0 < M K o n o < 

1.8 GeV/c 2 , which excludes the K*° mass region, yields S = -0.39lg;^(stat) and A = 
+0.10 ± 0.51(stat). The result on the S term is consistent with the previous result obtained 
for the decay B° -> fC 7 (K*° -> K° s tt°), S = -0.79lg;|g(stat) ± O.lO(syst), where 0.8 < 
Mk tt° < 1-0 GeV/c 2 was required. 

As discussed in Section I, to a good approximation, the SM predicts S = — £/sin20i for 
the B° — > KgKgK decay. Figure 3 summarizes the sin 20i determination based on Belle's 
S measurements using modes dominated by the b — > s transition [10]. For each mode, the 
first error shown in the figure is statistical and the second error is systematic. We obtain 
sin 20! = +0.39 ±0.11 as a weighted average, where the error includes both statistical and 
systematic errors. The result differs from the SM expectation by 2.7 standard deviations. 

IV. SUMMARY 

We have performed a new measurement of CP-violation parameters for B° — > K° S K®K° S 
decay based on a sample of 275 x 10 6 BB pairs. The decay is dominated by the b — > s flavor- 
changing neutral current and the KgKgKg final state is a CP eigenstate. Thus it is sensitive 
to a new CP-violating phase beyond the SM. The result differs from the SM expectation 
by 2.8 standard deviations. The combined result with the decays B° — > (f)K°, K + K~Kg, 
fo(980)Kg, rj'Kg, ooKg and i^^vr , which are also dominated by the b — > s transitions, differs 
from the SM expectation by 2.7 standard deviations. 
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We have also measured the time-dependent CP asymmetry in the decay B° — > i^^7r°7, 
which is also sensitive to physics beyond the SM. The invariant mass of the K^tt system is 
required to be between 0.6 and 1.8 GeV/c 2 , which is an extension of the previous analysis 
performed with the decay B° -> K*°-f (K*° -> K° s n°). The statistical error is much reduced 
from the previous analysis by including more events, and the result is consistent with the 
previous analysis. 

In both cases, measurements with a much larger data sample are required to conclusively 
establish the existence of a new CP-violating phase beyond the SM. 
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TABLE I: The event fractions wrong-tag fractions Wi, wrong-tag fraction differences Awi, and 
average effective tagging efficiencies e 1 ^ = e;(l — 2wi) 2 for each r interval for the DS-II. The errors 
for wi and Awi include both statistical and systematic uncertainties. The event fractions are 
obtained from J/ipK^ data. 



1 


r interval 


e/ 


Wl 


Awi 


e eff 


1 


0.000 


- 0.250 


0.397 ±0.015 


0.464 ± 0.007 


±0.010 ± 0.007 


0.002 ± 0.001 


2 


0.250 


- 0.500 


0.146 ±0.009 


0.321 ± 0.008 


-0.022 ±0.010 


0.019 ±0.002 


3 


0.500 


- 0.625 


0.108 ±0.008 


0.224 ±0.011 


±0.031 ±0.011 


0.033 ± 0.004 


4 


0.625 


- 0.750 


0.107 ±0.008 


0.157 ±0.010 


±0.002 ±0.011 


0.051 ± 0.005 


5 


0.750 


- 0.875 


0.098 ± 0.007 


0.109 ±0.009 


-0.028 ±0.011 


0.060 ± 0.005 


6 


0.875 


- 1.000 


0.144 ±0.009 


0.016 ±0.005 


±0.007 ± 0.007 


0.135 ±0.009 



TABLE II: The estimated signal purity and the signal yield N S i g in the signal region for B° — > 
K^K^K^ and K^^j decays. For the B° -> K%K%K% decay, results both before and after vertex 
reconstruction ( the latter in parentheses) are shown as events that do not have vertex information 
are also used to extract the direct CP violation parameter A. For the B° — > Kg7r 0, y decay, the 
result after vertex reconstruction is shown. 



Mode 


purity 


iVsig 




0.53 (0.54) 


88 ± 13 (63 ± 9) 




0.56 (0.56) 


72 ± 10 (54 ± 8) 




0.40 (0.40) 


16 ± 8 (8 ± 4) 



^vr u 7 0.50 105 ± 14 



TABLE III: Results of the fits to the At distributions. The first error is statistical and the second 
error is systematic. 

Mode SM expectation for S S A 

K° S K° S K° S -sin 2^ ±1.26 ± 0.68 ± 0.18 ±0.54 ± 0.34 ± 0.08 

i^7r°7 -2(m s /m 6 )sin 2<f>i -0.58^;^ ± 0.11 ±0.03 ± 0.34 ±0.11 



TABLE IV: Summary of the systematic errors on S and A. 





S K°K°K° 








Vertex reconstruction 


0.02 


0.05 


0.05 


0.06 


Flavor tagging 


0.04 


0.01 


0.01 


0.01 


Resolution function 


0.12 


0.05 


0.04 


0.04 


Physics parameter 


0.01 


0.01 


0.01 


0.01 


Possible fit bias 


0.03 


0.02 


0.02 


0.01 


Background fraction 


0.10 


0.06 


0.03 


0.04 


Background At shape 


0.08 


0.05 


0.01 


0.04 


Tag-side interference 


0.02 


0.01 


0.02 


0.06 


Total 


0.18 


0.11 


0.08 


0.11 



14 




30 




AE (GeV) AE(GeV) 



FIG. 1: The M bc distributions within the AE signal region for (a) B° -» K^K^K^ and (b) 
and the AE 1 distributions within the Mb c signal region for (c) B° — > iT^if^ifg and 
(d) -B —s- K^7r°7. Solid curves show the fit to signal plus background distributions, and dashed 
curves show the background contributions. 
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(b) i^7T°7 
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FIG. 2: The asymmetry, A, in each At bin with < r < 0.5 (top) and with 0.5 < r < 1.0 (bottom) 
for (a) B° — > A'^A'gA'g and (b) 5° — » Ag7r 7. The solid curves show the result of the unbinned 
maximum-likelihood fit. The dashed curves show the SM expectation with (S,A) = (— sin20i = 
-0.73, 0) for B° -» K° S K^K° S and with (S, A) = (0, 0) for B° -> A^7r°7. 
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Belle Results on Time-dependent CP Violation in b— ►s 




Average 

b— >s penguin 



sin2(|) 1 (WA) 



_l_ 



_J_ 



H 



H 



0.06±0.33±0.09 
0.49±0.18^ 
-0.47±0.41±0.08 
0.65±0.18±0.04 

0.75±0.64t°^ 

0.30±0.59±0.11 

-1.26±0.68±0.18 

0.39±0.11 

0.726±0.037 
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FIG. 3: Summary of sin2</>i measurements performed with B decay governed by the b — > sqq 
transition. The world-average sin2</>i value obtained from B — ► J/tpK° and other related decay 
modes governed by the b — » ccs transition [7] is also shown as the SM reference. 
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